Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN I^TeX style file v2.2) 



Gravitational wave background from coalescing compact 
stars in eccentric orbits 

V.B. Ignatiev\ A.G.Kuranov^, K.A. Postnov^'^, M.E. Prokhorov^ 

^ Faculty of Physics, Moscow State University, 119899 Moscow, Russia 

^ Sternberg Astronomical Institute, 119899 Moscow, Russia 

^ Max-Planck Institut fiir Astrophysik, 85740 Garching, Germany 



Accepted 2001... Received 2001 ... 



ABSTRACT 

Stochastic gravitational wave background produced by a stationary coalescing 
population of binary neutron stars in the Galaxy is calculated. This background is 
found to constitute a confusion limit within the LISA frequency band up to a limiting 
frequency vum ~ 10^'^ Hz, leaving the frequency window ^ lO^'^-lO^^ Hz open for 
the potential detection of cosmological stochastic gravitational waves and of signals 
involving massive black holes out to cosmological distances. 
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1 INTRODUCTION 

In a short time, gravitational wave astronomy is expected 

■ to start collecting first data from different astrophysical 
\ sources of gravitational radiation, most probably from co- 

■ alescing binary neutron stars and black holes (for a recent 
review, see Grishchuk et al. 2001). The first large gravita- 
tional wave (GW) interferometers LIGO, VIRGO, GEO-600 
and TAMA-300 will be sensitive to dimensionless metric 
variations in the frequency range 10-1000 Hz at a level of 
h ~ lO"'^^ (Braginsky 2000 and references therein). Due to 
seismic noise limitations lower frequencies (10~^ — 10"^ Hz) 
can be studied only with space-born antennas, such as LISA 
laser interferometer (Bender et al. 2000) . Prospects for find- 
ing astrophysical sources within the LISA frequency band 
are very good, especially for supermassive binary black holes 
which can be observed by LISA with a record signal-to noise 
ratio of order 1000 from cosmological distances (Vecchio et 
al. 1997). 

In addition to signals from individual sources, GW 
interferometers can detect stochastic gravitational waves. 
Such waves (GW backgrounds, or GW noises) can be pro- 
duced by a large collection of unresolved individual sources, 
e.g. by binary stars (Mironovskij 1965, Lipunov and Postnov 
1987, Bender and Hils 1997, Kosenko and Postnov 1998) or 
rapidly rotating neutron stars (Giazotto et al. 1997, Post- 
nov and Prokhorov 1997, Ferrari et al. 1999). These GW 
backgrounds are usually considered as unwanted additional 
noises to the intrinsic noises of GW detectors, since they 
could potentially mimic stochastic GW backgrounds of cos- 
mological origin (primordial or relic GW) that bring a valu- 



able information about physical processes in the very early 
Universe (Grishchuk et al. 2001 and references therein). 

Ordinary galactic binaries, in which at least one of 
the component is a normal main-sequence star, mostly 
contribute at low-frequency band {v ~ 10~^ — 10~^ Hz) 
(Mironovskij 1965). At higher frequencies 10"* — 10~^ 
Hz compact coalescing white dwarf (WD) binaries should 
form a noticeable background above the LISA noise curve 
(Lipunov and Postnov 1987, Bender and Hils 1997, Gr- 
ishchuk et al. 2001). Generally, the level of the GW noise 
within a fixed frequency interval from a collection of N un- 
resolved independent sources in terms of dimensionless am- 
plitude is /i oc %/iV. In case of binaries that loose the orbital 
angular momentum due to gravitational radiation back re- 
action, this level is /i oc y/lZ, where TZ is the coalescing rate 
of these binaries (see Grishchuk et al. 2001 for more detail). 

Another important quantity is the upper frequency vum 
above which each individual source can be resolved dur- 
ing one-year observation time (i.e. within the frequency bin 
Ai^ = 3 X 10~** Hz). In fact, it is this limiting frequency 
which mostly relate to the challenging task of detection of a 
relic GW background: at > we will be able in princi- 
ple to single out individual sources and thus have prospects 
to register cosmological GW noise using one interferometer 
(Grishchuk et al. 2001). 

The obvious condition for vum reads TZAT{i'iim) — 1, 
where AT(i^) is the time each source 'spends' within the 
frequency bin at a given frequency ly. For example, in the 
simplest case of a collection of galactic binary WD in cir- 
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cular orbits which coalesce at a constant rate TZ = 7?.3oo x 
(1/300) yr-i 

/^i>™(WD) ^ (1.2 X 10-^Hz) n%f 




where M is the chirp mass of the binary. For merging binary 
neutron stars (NS) vum calculated as above would give even 
smaller value, i'Hm(NS) = 3 x 10~* Hz, since the binary NS 
coalescence rate in most optimistic scenarios isTZns ~ 10~* 
yr"^ (Grishchuk et al. 2001, Kalogera et al. 2000). Note that 
the uncertainty in this rate by a factor of 3 or so is of minor 
importance since viim oc TZ^^^^ . 

However, there is an important difference between the 
merging WD and NS. Binary WD must have almost circu- 
lar orbits from the very beginning since they result from 
a spiral-in process during the common envelope stage. Un- 
like binary WD, binary NS must have (and this is what we 
actually observe in the known binary radiopulsars with NS 
companions) a significant eccentricity at birth, since they are 
formed after two supernova explosions with substantial mass 
loss from the system. The possible asymmetry of supernova 
explosion leading to the natal kick velocity of newborn NS 
additionally affects the orbital parameters leading to higher 
orbital eccentricities. 

It is well known (Peters and Mathews 1963) that an ec- 
centric binary system emit GW in a wide frequency band. 
This means that the high-order harmonics from an eccentric 
binary system should be observed at frequencies u > 2i/k 
{2vk is twice the Keplerian orbital frequency at which a cir- 
cular binary radiates GW), so the entire population of galac- 
tic eccentric binary NS should contribute in a broader fre- 
quency band than analogical circular binaries would do. The 
effect may be not small, since the total number of NS-I-NS 
binaries in the Galaxy is quite large, as the simple estimate 
below shows. Assume the stationary formation rate TZf oi 
NS binaries in the Galaxy (which is a reasonable approx- 
imation during at least the last 5 billion years). Then the 
total number of them is of order TZf x Tgai, where Tgai is 
the galactic age. Note that the actual NS formation rate is 
higher than the NS merging rate TZns since a NS binary 
could merge only if its initial semi-major axis at and ec- 
centricity ei satisfy the condition Tcoai{ai,ei) < tgai, where 
Tcoai is the time the systems takes from the formation to 
coalescence due to GW orbit decay. So the minimum num- 
ber of NS binaries in the Galaxy is estimated to be of order 
TZns x Tgai = 10"* x 10^° ~ 10^ 

To calculate the GW background produced by binary 
NS, the following steps should be done. (1) NS binaries form 
with some initial distribution over orbital semi-major axes 
and eccentricities Fin{a,e). This distribution can be found 
from evolutionary calculations using e.g. binary population 
synthesis method (Lipunov et al. 1996). (2) Next, if orbital 
parameters of such binaries evolve only due to GW emission, 
a stationary distribution function Fst{a,e) can be recovered 
(Buitrago et al. 1994). (3) Knowing the GW spectrum of 
a binary with given (a,e) (Peters and Mathews 1963), the 
final GW noise from these sources can be computed. 

In the present paper we specially focus on the effect 
of eccentricity of coalescing NS binaries on the level and 
limiting frequency of the stochastic GW background they 
form. We find, as expeted, that the level of this background 




Figure 1. The initial binary NS distribution Fi„(a,e) for the 
kick velocity amplitude Vf; = 200 km/s. Initial semimajor axis a 
is in solar units. 



is lower than from coalescing galactic binary WD, and the 
limiting frequency i^umie 7^ 0) is higher than calculated for 
circular case. The initial (model-dependent) distribution of 
NS binaries only slightly affects the result. This leaves the 
frequency interval ~ 10~^ — 10~^ Hz open for searches of 
primordial GW backgrounds. 



2 INITIAL GALACTIC BINARY NS 
DISTRIBUTION 

To calculate the initial binary NS distribution Fin{a,e) 
in the Galaxy we used the population synthesis method 
described by many authors (Lipunov, Postnov, Prokhorov 
1996; Lipunov, Postnov, Prokhorov 1997; Portegies Zwart, 
Yungelson 1998). It uses the simulation of evolution of a 
large number of binaries with initial parameters (masses 
of the components, semimajor axes, eccentricities etc.) dis- 
tributed according to some (taken from observations or 
model) laws. There are also evolutionary parameters, such 
as efficiency of common envelope stage, kick velocity dur- 
ing supernova explosion, initial spin of compact objects, etc. 
(see Grishchuk 2001 for more detail), of which kick velocity 
imparted to compact object (neutron star or black hole) 
at birth mostly impacts the resulting distribution Fi„{a,e). 
We assumed a Maxwellian form of the kick velocity distri- 
bution and varied its amplitude Vk from to 400 km/s. 
The resulting normalized initial distribution Fi„{a,e) for 
Vk = 200 km/s is shown in Fig. ^ 

NS binaries are formed in a very broad interval of a 
and e, but interesting for us here will be only those that 
can coalesce over the galactic age Tcoai < Tgai, since only 
such binaries can form a stationary distribution. It is seen 
from Fig. |l| that most binaries form with small or moderate 
eccentricities (e < 0.5) in a wide range of semimajor axes 
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Figure 2. The stationary distribution function Fst for binary NS 
in eccentric orbits which coalesce due to GW emission in a time 
interval shorter than the galactic age 10^'^ years. 



Clin ~ 1 — 3OOi?0. Increase in the kick velocity amplitude 
broadens the initial distribution function over eccentricity. 

3 STATIONARY BINARY NS DISTRIBUTION 

GW back reaction causes orbital decay of binary system. For 
an eccentric binary star, this can be described in quadrapole 
approximation by ordinary differential equations for orbital 
semimajor axis and eccentricity (Peters & Mathews, 1963): 
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Given the initial function Fin{a, e), it is straightforward 
to find a stationary distribution function (a, e) (Buitrago 
et al. 1994): 
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At present, only some fraction of systems from the ini- 
tial distribution with Tcoai < Tgai reaches stationarity. The 
function FatiyKjC) is calculated for this part of Fin{a,e) 
from Fig. |l| and is shown in Fig. | 2| In this figure we use the 
Keplerian frequency uk = \/ (l/27r)G(mi + m2)/a^ instead 
of the semi-major axis a. 



4 GW SPECTRUM FROM NON-CIRCULAR 
BINARY STARS 

In the simplest case of two point masses in a circular orbit, 
the energy losses caused by the quadrupole gravitational 



wave emission Lq = ^^('^^) defined by (Einstein, 1916 
1918; Landau & Lifshiz 1971): 

32 G* m?mi(mi -f 7712) 



Lo = 
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(5) 



where G is Newton's gravitational constant, c is the speed of 
the light in a vacuum and m\,m2 are the masses of stars. All 
GW radiation in this case is emitted at a single frequency 

For a non-circular orbit with eccentricity e the GW lu- 
minosity increases (Peters & Mathews 1963): 
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The emission now is widely spread over higher-order har- 
monics to the main frequency (Peters & Mathews 1963): 
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Here J„ is Bessel function and n is the number of the har- 
monics. 

The total energy emitted per second in a frequency in- 
terval Al/ by all binaries in the Galaxy is 
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where Fst(^,e) = Fst{vk,e) is a stationary distribution 
function of binary NS. In the last equation we take into 
account Av <^ v. The amplitude of high-order harmonics 
rapidly decreases, so we stopped the summation of harmon- 
ics for n > niim, num ■ g{niim,e) = e x max„[g(n, e)]. We 
assumed e — 10~*, 10~®, lO"'^". Decreasing e increases the 
total number of harmonics which contribute to the given 
frequency bin, but practically does not change the number 
of the most powerful harmonics within it (see Fig. 7). 



5 STATIONARY STOCHASTIC GW NOISE 
FROM COALESCING BINARY NS 

First, consider the GW spectrum from a delta-function ini- 
tial distribution: Fin{a,e) = S{a — ao,e — eo), that is assume 
all binaries to form at one point (ao, eo) in (a — e) configu- 
ration space at a rate of TZns = lO"** yr^^. The spectrum 
formed by the first three harmonics and the total station- 
ary GW background for this case are shown in Fig. ^ and 
Fig. ^for ao — 5Rq and eo = 0.5, 0.9, respectively. It is seen 
that increase in the initial eccentricity strongly affects the 
shape of the spectrum up to some frequency at which the ec- 
centricities of binaries in the stationary distribution become 
sufficiently small. Above this frequency only the second har- 
monics from almost circular binaries contributes to the total 
spectrum. The non-monotonic dependence of the spectrum 
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Figure 3. The resulting GW spectrum and tlie contribution from 
the first, second and third harmonics. All systems form with ao = 
5-R0, eo = 0.9. For comparison we show the spectrum for circual 
NS binaries coalescing with the same rate, and with initial orbital 
periods longer than 2 X lO^s. 



1 ^ yr 1 




log i/.llz 



Figure 4. The same as in Fig|3|for systems with ao = 5-Rq, eo ■ 
0.5. 



formed by the harmonics is due to the non-monotonic be- 
haviour of the energy emitted at each harmonics with ec- 
centricity (Eq. (|)-(|)). 

Now we turn to the expected GW stochastic back- 
ground from galactic NS binaries. We shall assume all the 
sources to lie at one distance r = 7.9 kpc, which is close 
to the average distance to stars in our Galaxy. This simpli- 
fying assumption is unlikely to change our general conclu- 
sions. At each frequency v we sum up the GW flux from 
all the harmonics that fall within the chosen frequency bin 
Afc' = 3 X 10~* Hz from lower-frequency non-circular systems 
in the calculated stationary distribution Fst{v,e) (Eq. (^)), 
as modified by including only the part that has come to equi- 
Ubrium in the lifetime of the galaxy, as discussed earlier. 

The resulting noise curve is shown in Fig. |5| in terms of 
dimensionless amplitude h 



Figure 5. GW background from stationary binary NS population 
calculated for a model spiral galaxy with the NS-I-NS coalescing 
rate T^jvs = 10~* yr~^. The limiting frequecy above which the 
most powerful harmonics can be resolved in one- year observation 



IS Uli 



0.8 mHz. If no harmonics were subtracted from the 



frequency bin 3 X 10~° Hz, the background would continue to- 
ward higher frequencies as shown by the long-dashed line. THe 
lower dashed line shows contribution from extragalactic NS-(-NS 
binaries. The upper solid line schematically represents WD-I-WD 
confusion limit from Bender and Hils (1997) iyu^ ~ 1.2 mHz, 
Eq. (1)). The proposed LISA sensitivity is from Thorne (1995) 
[Fig. 14] 
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As expected, the NS-f-NS confusion noise lies below 
WD-fWD curve, mainly due to lower TZns- High-order har- 
monics from non-circular NS binaries mostly contribute at 
lower frequencies, so starting from u ~ 10"* Hz the calcu- 
lated noise curve practically coincides with that formed by 
circular NS binaries coalescing with the same rate 7?.jvs = 
10-^ yr-i. 

More important in the non-circular case is the increase 
of the limiting frequency uum above which individual sources 
can be resolved during a one-year observation. To estimate 
it, we calculate the number of harmonics inside the fre- 
quency bin Au. Formally, this number is more than one at 
any frequency, but at high frequencies the main contribution 
must come from circular binaries (2-nd harmonics), the total 
GW power from higher-order harmonics from low-frequency 
eccentric systems becoming gradually less and less. So the 
number of harmonics contributing in a given frequency bin 
Au = 3 X 10~* Hz were counted, starting with the strongest 
one and continuing until 99bin (Eq.(^) had been included. 
The number of such harmonics as a function of frequency 
and the assumed kick velocity amplitude N{iy, Av) is shown 
in Fig. ^. The limiting frequency vum is determined from 
equation N{i'iim, Au) = 1. To within uncertainties of our 
calculations {TLns, Fi„{a,e), etc.), vum — 10~^ Hz, close 
to the break in the confusion limit from binary WD. For 
comparison, in the same Fig. ^ we show N{u, Au) for sta- 
tionary circular NS binaries coalescing with the same rate 
7^Jvs = lO""* yr"^ {ua^ie = 0) cx 3 x 10"'' Hz). The effect 
of increasing the chosen level for the GW power from the 
strongest harmonics in the bin is illustrated in Fig. ^. In- 
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Figure 6. Number of harmonics from non-circular galactic NS 
binaries producing the stochastic noise in the frequency bin 
Au = 3 X 10~* Hz as a function of the kick velocity. The curves 
are labeled with kick velocity amplitudes 100, 200 and 300 km/s 
assumed in the population synthesis calculations. The straight 
line is for the circular NS binaries. The limiting frequency un^y^ 
above which individual harmonics can be resolved follows from 
the condition N{uii^, Au) = 1. 




Figure 7. Number of the most powerful harmonics from non- 
circular galactic NS binaries which contribute 95%, 99% and 100% 
of the total GW energy collected in the bin Ai/ = 3x 10~* Hz. The 
effect of increase in the number of harmonics njim determining 
the total energy collected in the bin (Eq. (9)) is shown for the 
total number of harmonics (100%) for different e. Kick velocity 
amplitude 200 km/s. 



creasing the level from 95% to 100% changes the limiting 
frequency vum by almost an order of magnitude. 

This procedure, however, is not complete. Indeed, con- 
sider a frequency bin next to thus defined viim- We may ask 
the question which harmonics will be most probably absent 
inside it. The answer is those which are the least probable 
at this frequency. The probability to find a harmonics in 
the bin at a given frequency is determined by the number 
of the harmonics and the stationary distribution function 
of sources /(i/, e) It happens that this is harmonics number 



i 



one (mostly due to a steep decrease of the systems' station- 
ary distribution function) . If we remove all 1st harmonics of 
systems with orbital frequencies falling into the chosen bin 
from the sum Eq.(^), we are left again with some (smaller) 
GW power in the bin and may wish to find the limiting fre- 
quency exactly in the same way as described above (i.e. by 
fixing some level and summing up the strongest harmonics 
up to this level), call it viim,-i- Above this new limiting fre- 
quency we can repeat the entire procedure to find uum.-2 
(this time the 2d harmonics is least probable to be found in 
the bin next to the new limiting frequency), etc. until the 
noise level of the detector is reached (it is less interesting for 
us to see the changes in the spectrum below it) . The step-like 
fine which continues the spectrum above ~ 1 mHz in Fig. 

illustrates this procedure. This line of course do not rep- 
resent the "real" GW noise curve from binary NS and just 
give a feeling of how it most probably behaves at uum, above 
which individual harmonics from coalescing binary NS can 
be singled out. 

Extragalactic NS binaries would also form an isotropic 
confusion noise, but the level of extragalactic GW back- 
grounds is generally an order of magnitude lower than the 
galactic one (e.g. calculations of Kosenko and Postnov, 1998) 
and is beyond reach by the expected LISA sensitivity (the 
bottom dashed curve in Fig. 5). The limiting frequency uum 
for extragalactic NS binaries finds from the condition for cir- 
cular systems and can be as higher as ~ 0.2 Hz (Ungarelli 
and Vecchio, 2000). 



6 CONCLUSIONS 

We studied the effect of non-circularity of realistic coalesc- 
ing NS binaries on stochastic galactic background formed 
by stationary NS-I-NS distribution Fst{i', e). The level of the 
background is found to be by one order of magnitude less 
than from coalescing binary WD systems, in correspondence 
with the lower NS binary coalescence rate. The limiting fre- 
quency above which individual harmonics can be resolved 



at a 99% level over this noise is vir 



0.8 X lO-" Hz, 



close to viim for coalescing binary WD. Detection of an 
isotropic stochastic GW signal by LISA at higher frequency 
would strongly indicate its cosmological origin. So our study 
confirms that within the frequency range 10-=*-10-^ Hz 
there are prospects to detect cosmological stochastic GW 
by means of one space interferometer LISA as suggested in 
Grishchuk et al 2001. 
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